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Abstract 

We study the quantum coherent-tunneling between two Bose-Einstein con- 
densates separated through an oscillating trap potential. The cases of slowly 
and rapidly varying in time trap potential are considered. In the case of 
a slowly varying trap we study the nonlinear resonances and chaos in the 
oscillations of the relative atomic population. Using the Melnikov function 
approach, we find the conditions for chaotic macroscopic quantum-tunneling 
phenomena to exists. Criteria for the onset of chaos are also given. We find 
the values of frequency and modulation amplitude which lead to chaos on 
oscillations in the relative population, for any given damping and the non- 
linear atomic interaction. In the case of a rapidly varying trap we use the 
multiscale expansion method in the parameter e = 1/^2, where i7 is the fre- 
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quency of modulations and we derive the averaged system of equations for 
the modes. The analysis of this system shows that new macroscopic quantum 
self trapping regions, in comparison with the constant trap case, exist. 
PACS numbers: 03.75.Fi, 05.30.Jp 
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I. INTRODUCTION 



It has recently been shown that there exists a macroscopic quantum phase difference in 
processes connected with atomic waves. Namely, this effect is observed between two tunnel- 
coupled Bose-Einstein condensates [Q- 0. The coupling occurs due to a trap potential 
which has the form of a double- well potential. The theoretical investigation shows that here is 
possibly an interesting phenomenon of periodic oscillations of the atomic population between 
condensates and quantum self-trapping of population in the dependence of the relative phase 
between condensates Analogous phenomena which have been studied are the ac 

Josephson effect ^ and the periodic exchange of power and switching of electromagnetic 



waves between cores in nonlinear optical couplers |10 



In particular, there is a direct analogy between the tunneling phenomena in two prolon- 
gated Bose-Einstein condensates and two tunnel-coupled single mode optical fibers. In this 
optical analogy, the role of the chemical potential is played the propagation constant, and 
to the nonlinear interaction between atoms is analogous to the Kerr nonlinearity of optical 
media. The tunnnel-coupling arising from the overlaps of the electromagnetic fields outside 
of dielectric cilinders (fibers) |[TT|-|r3| exactly corresponds to the tunnel coupling between 
two Bose-Einstein condensates due to the overlaps of the wavefunctions. 

In this context , it is natural to investigate the influence of a time varying coupling 
on the quantum coherent tunneling process. In the nonlinear optical coupler analogy this 
corresponds to the variation in the longitudinal direction coupling As was shown 

recently, the variation in time of the trap potential can lead to resonant oscillations of the 



Bose-Einstein condensate The numerical solution of the Gross-Pitaevsky equation 

shows that the rapid variation in time of the trap potential can lead to a bifurcation of the 
effective (averaged ) form of the trap potential and consequently to splitting of the BEG 

From what we have said above, we could wait for new phenomena in the process of quan- 
tum coherent atomic tunneling (QGAT) process too. Recalling the analogy with Josephson 
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effects in coupled junctions, we note that in the last case it is very difficult to produce such 
kind of effects, the difficulty coming from the fact that it is hard to implement the time 
varying overlap properties of junctions. 

Thus, in this article we will study the atomic tunneling between two tunnel-coupled 
EEC's in a double-well time dependent trap. Namely, in Section II we will study the 
nonlinear resonance phenomena in the interference and, in particular, the stationary regimes, 
with damping taked into account. The phase and population damping are considered as 
well. The phase locked states are interesting because the suppression in these states of 
the fluctuations of the relative phase between two condensates. The regions of regular and 
chaotic oscillations of the relative atomic population are calculated in Section III, using the 
Melnikov function approach. In Section IV the macroscopic quantum interference in two 
overlapped condensates with rapidly varying traps is investigated. We derive the system of 
averaged equations for coupled modes and define in Section V new regimes of parameters 
for the macroscopic quamtum self-trapping (MQST). 

II. FORMULATION OF PROBLEM. SLOWLY VARYING IN TIME TRAP 

The problem of Bose- Einstein Condensates in a two double- well time dependent trap can 
be described by the following model 

th^ = {Ei{t) + ail^'iH^'i - K{t)tlj2 

ih^ = {E,{t) + a2^2n^2 - K{t)^,, (1) 

where the parameters Ei,ai, K{t) are defined by the overlaps integrals of the time dependent 
Gross- P it aevsky eigenfunctions as given in [Q. are parameters of nonlinear interactions 
between atoms and ctj ~ (70 = 4:TTh'^a/m, a is the atomic scattering length. The attractive 
case can be obtained analogously, using symmetries of the equation This system of equations 
is valid in the approximation of a weak link between condensates. Comparison with the 
numerical solution of GPE shows that the system (|l]) is a reasonable approximation for 
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z < 0.5 — 0.6. For the strongly overlapped condensates the system should be modified. For 
a periodically varying K{t) = Kq + Ki sin(nt) we have limitations on the parameters Ki 
and VL for which the two-mode approximation is valid. The small deformation of bound 
states requires Ki <^ Kq. To avoid the resonances of the modulations with normal modes of 
the spherical well we should require that the difference of energy between the ground state 
and the first normal mode be larger than the modulation energy, i.e. 5E ^ hVt. Taking 
into account that 5E ~ Tiuoq, 0, where uoq is the harmonic oscillator frequency, we find the 
condition <^ uoq. Taking into account that cuq ^ = Sii'o ^ the frequency of the linear 
oscillations of the atomic population, defined by the tunneling frequency, we conclude that 
by a proper choice of parameters we can have the resonant case f2 ^ cuj;, or the rapidly 
varying modulation uo^ . 

Introducing new variables Ui = y/Wiexp^iOi), z = {Ni — N2) /Nt, Nt = A^i + N2^'ip = 
61 —62, where Ni, and 6i are respectively the number of atoms and phases in i-th trap, we 
get the following system 



zt = -2K{t)Vl - z^sin^ - r/$t, (2) 
2K(t)z 

^t = JyAz+ , ^ \ cos$ + AE(t), (3) 

where A = {ai + a2)NT, and z/ = ±1 for the positive and negative atomic scattering length 
respectively. In what follows, we will consider the case z/ = 1, (a^ > 0), if not stated 
otherwise. We have included in Eq(0) the damping term //^^(t). It appears if we take 
into account a noncoherent dissipative current of normal-state atoms, proportional to the 
chemical potential difference Afi. For the other type of overlapping condensate there might 
be another type of damping. For example, for the two interacting condensates with different 
hyperfine levels in a single harmonic trap the damping has the form i]z{t). 

The Hamiltonian of the unperturbed system {K{t) = const, Ai? = const, 77 = 0) is 



Az 



2 



H = — - 2K^l - ^2 cos(<l') + AEz . (4) 

Let us consider the case of periodic modulations of the tunnel coupling parameter, when 
K{t) = Kq + Ki sin Qt and the time-periodic energy difference is AE{t) = AE-\-AEi sin{Qt). 



It is useful to introduce the parameter A = $ — Qt, which is a slow-varying function of the 
time in a comparison with the period Q. The averaged system of equations can be found by 
averaging over the period of fast oscillations 27r/r2. Doing so, we get the system 

Zt = —Ki cos AVl — z'^ — i]At — 1]Q, (5) 

KtZ - 

At = -n + AE + Az . _ sin A. (6) 

Let us now look for the fixed points of Eq.(^. Introduce first the notation a = f2 — AE. 
From the system (1^) we find the following equation for the fixed points 



1 (7) 

sin 2 Ac, 



For small Zc, z"^ << 1 we can find the explicit solution 

n-AE 



Ac = arccosf 



(9) 

A - ^Kf - r]^D? 

The existence of a stationary solution is connected with the circumstance that the damp- 
ing of the oscillation in quantum tunneling is compensated by the periodic variation of the 
trap. In reference the influence of the damping on the macroscopic quantum interference 
has been investigated numerically and degradation of the interference due to damping has 
been analyzed. The predicted stationary solutions, where these effects are compensated 
by the variable driving, are indeed phase - locked solutions. Thus we can think that the 
fluctuations of the relative phase will be suppressed in such states. It is an important predic- 
tion since EEC's have phase fluctuations and it is difficult to maintain the constant relative 
phase of condensates. This phenomenon can be useful to develop the phase standard for the 
Bose- Einstein condensate- a problem attracting considerable attention recently 

We can estimate the value of fixed points for an experimental situation. For example, 
when Vl = 0.7,1] = 0.1, Ki = 0.2, A = 2.5, we obtain the fixed points zi^ = 0.2599, 2:20 = 
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0.3037. Stability analysis shows that only the first point is stable and corresponds to the 
global attractor. The results of the numerical simulations of the initial system for the 
relative population z{t) are presented in Fig.l. It shows that the agreement between theory 
and numerical simulations is good. 

The same system of equations also describes the nonlinear Josephson type oscillations 
in the relative population of a driven , two-component Bose-Einstein condensate. The two 
component condensates where two different hyperfine states can be populated and confined in 
the same trap. A weak driving field couples two internal states and as result the oscillations 
in the relative population will occur In this case the damping term has the form 

r = —riz. Repeating the above mentioned procedure we obtain the averaged system 



Zt = —Ki \/l — z"^ cos A — r]z, 

K z - 

At = Az + AE-n ^ ^ _ sin A. (10) 



z 



In the model |20|, AE corrsponds to /i2 — /ii — f^, where fii are chemical potentials and k 
is frequency detuning k = cud — ujq, uq is the separation in frequency of two hyperfine states, 
Uq is the frequency of drive. The fixed points are: 



arccos 
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Az -a S^ri ^ rill 



Here a = Vt — AE. For the case of small 2;^ << 1 we have the estimate for Zc 

AE-n , {AE-ny{Kf-r^') 
A-K, 2Kr{A-K,Y ' 

The typical values of parameters aie:Q = 0.7, /i2 — — n = AE = 1,1] = 0.1, i^i = 
0.2, A = 4, and the fixed points are z^c = —0.111, z^c = —0.13. The analysis shows that the 
first fixed point is stable and second is unstable. 

In Fig. 2. we present the plot of solution of the full system for the damping term r]z with 
the above mentioned set of parameters. All solutions decay to the first fixed point. 



III. CHAOTIC OSCILLATIONS IN RELATIVE ATOMIC POPULATION 



In this section we consider the interference in the system which in some hmit is mathe- 
matically equivalent to the motion of a particle in a double-well potential under parametric 
periodic perturbations. The particle' motion is of course unharmonic and the nonlinear res- 
onances between oscillations of a particle and oscillations of the trap potential are possible. 
One of consequences of this is the possible appearence of chaotic dynamics in such systems. 
We will now study the conditions for the appearence of the chaotic dynamics in atomic 
tunneling. 

To proceed it will be useful first to consider some aspects of the unperturbed system, 
that is, when the coupling is constant. In this case the system (H) is equivalent to the quartic 
(Duffing) oscillator []6|,|23|. 

We can write the equation for z{t) as 

— f . (13) 

where the potential V{z) is given by 

Viz) =z\a + bz'), a = 2Kl - _, 6 = _. (14) 

and H is the Hamiltonian. The total energy of the effective particle is 

z"^ 

Eo = ^ + V{z) = 2K^ - H^/2 . 

When E'o > we have the oscillating regime with < z{t) >= (nonlinear Rabi oscillations). 
This case corresponds to the periodic flux of atoms from one BEG to other. When Eq < 
the motion of effective particle is confined in the one of the wells. This case corresponds to 
the localization of atomic population in one of the condensates - the so called macroscopic 
quantum self-trapping (MQST). The case Eq = 0, i.e. H = 2K, corresponds to the separa- 
trix solution separating these two regimes. It is interesting to investigate the dynamics of 
atomic population, when the inital parameters ^^(0), $(0) are close to the separatrix of the 
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unperturbed systsem. From the general theory of nonhnear driven oscillations we can 
expect the appearence of chaotic macroscopic tunneling phenomena. 
The separatrix solutions for the right hand side well(condensate) are: 

Zs{t) = -i/ ^sech(\/2at) , 
V b 



. 9,^ , a^sech^(v^t) tanh^(v^t) 

sin^(<l>s) = — ^= , (15) 



where H = H, = 2K. 

The motion near the separatrix is very sensitive to the change of the initial condition and 
a stochastic layer appeas, leading to the chaotic dynamics of the solutions of the system @. 



The Melnikov function method is effective to find the regions of chaotic behavior ||21| , |22 
The system (Q) may written in the form = f + eg and has a hyperbolic fixed point. 
According to this method, we need to calculate the Melnikov function M(to) 

/oo 
{f2s9ls{t,to) - fls92s{t,to))dt. (16) 
-00 

As is well known, the existence of zeros of M{to)) indicates the intersection of separatrices 
and the existence of homoclinical chaos. Substituting the expressions for fi,gi from the 
system (^ and calculating the integrals, we find the final expressions for the Melnikov 
function. We shall give the results for two possible damping terms. 

A. The case of the damping in the form 7]^t 

The integrations in (p!6[) are cumbersome, and will be omitted. The final result for the 
Melnikov function is: 

Mi(to) = -i^iFiCOs(fito)-i^2r/ , (17) 

where 

'"4^o&sinh(^)' 
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where 

arctan(a / \J ah ~ a?' 



Gi = [-2h' {2Hs - A) {hHs - 4aA + 36A)] - . 

yah ~ a? 

Aa^K^ -bahk^ + 2h'^{AHl -AH.K + ^k^)] . (19) 
The simple zero of M(to) is absent if the condition 

V > T^, (20) 

is satisfied. In Fig. 3 we plot this criterion in the {Ki,Q) plane for rj = 0.1 . Regions below 
the curve correspond to regular oscillations and above to chaotic oscillations of z{t). The 
criterion gives a good lower bound on the region of chaos in this plane. We verified this 
lower bound in {Ki,Q) plane by numerical simulations and found a good agreement with 
the formulae (^Ol). 

When the damping coefficient is zero we have the expression for the width of stochastic 
layer near the separatrix as 
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D = „ ^. (21) 

4ifo&sinh(:'^^ ^ ^ ^ 



In Fig. 4 we plot the width of the stochastic layer as a function of the frequency for 
Ki = 0.2, A = 9.9, ^(0) = 0.6, $(0) = 0. The maximum of the width is given by the solution 
of equation 

, , vrfi , 4^ 

iZtanh — = . 

2V2a vr 

For this choice of parameters we have = 1-89. For > 2,a > 1 we get the estimate 
Qm ~ 4^/2a/7r. It can be seen that, for the high frequencies >> y/a, the stochastic layer 
is exponentially narrow and motion is regular. The analysis of this case will be performed 
in detail in the next section. In Fig.5 we plot the typical oscillations in the relative atomic 
population for A = 9.9, Kq = 1, Ki = 0.2, z{0) = 0.6, $(0) = for = 1.5. It turns out that 
for Q between 1.5 and 3.22 we have chaotic oscillations, and for small or large fl, regular 
oscillations. In Fig. 6 the phase portrait is given. In Fig. 7 the influence of the damping on 
chaos is illustrated. Parameters are the same as in Fig.5 and rj = 0.15. In this case the 
oscillations become regular. 
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B. The case of the damping term in the form riz{t). 



The calculation of the integrals which give the Melnikov function is again cumbersome, 
and one find the following result: 

M(to) = Fi cos{nu) + (22) 

where 



^ arctan f ) + ^A^. (23) 



!2{h-a) \^b-aj 2V2b 

The zeros of M{to) are absent when the condition 

V > f , (24) 

is satisfied. 

If we consider for example the case when Ki = 0.2, Q = 0.7, A = 2.5, z{0) = 0.6, $(0) = 
TT, the damping constant should he rj = 0.15. In Figs. (7,8) the chaotic and periodic oscilla- 
tions of the relative population in the dependence of the frequency for fixed -2(0), $(0) are 
plotted. 

IV. RAPIDLY VARYING TRAP POTENTIAL 

In the rapidly varying case K = K{^) we can apply he averaging techniques based on 
the multiscale expansions |]2^p5| . Let us search the solution in the form of the sum of the 



slowly varying components U, V and rapidly varying corrections Ui, Vi 

Mt) = U{tk) + eui(C) + 6^(0 + ... 

Mt) = V{tk) + ev2{C) + e^V2{0 + ... (25) 

where C, = t/e,e « 1,, tk = e^t are slow times. The averages over a period of rapid 
oscillations are < ui >=< Vi >= 0. The functions Ui,Vi are assumed to be functions of the 
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rapid time C and the slow varying functions U, V. The derivatives can be written in the 
form 

dU dU dU 

— = — + €— + ... (26 

dt dto dti 

and analogously for the derivative of V. 

Below, we will assume Ei — E2,. Then the time dependence in E(t) can be removed by 
a simple field transformation u{v) = u{v) exp{{i/h) J E{t')dt'). In the zeroth order in e we 
get the equations 

+ = ^«o-i^(CH + a|«o|'«o (27) 
oiq ac, 

^^(^ + ^) = - ^(0«o + a\vo\\o (28) 
Requiring the exclusion of the secular in ^ terms, we find the corrections Ui, Vi 

Ui = j^voini- < Hi >) (29) 
vi = ■^Mo(a*i- < A*i >) (30) 
where fii = Jq{K{s)— < K >)ds and the system for uq^vq becomes: 

ih—^ — Euq— < K > vq + aluol^Uo (31) 
dvo 

ih—-^ — Evq— < K > uo + alvol'^Vo (32) 
In the order e we have 

and we obtain, for the second order corrections, U2,V2, the elimination of the secular terms, 
results in: 

i ^rdVo i r I i2 2 *m 

-[{fxi- < fxi >f - 2M]uo , (34) 



-[(/Xi-</Xi>)2-2MH , (35) 
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where 

fJ'2 — [ (/^iC~ < /^i >)ds, M = {< ii\> — < Hi > 
Jo 

^Prom this order we get the equations: 

[2\vo\\o-vlul] , (36) 



— IVId^I i/.^ — 7-^-. "I 

Rn^ O/irvA/T 

[2|Mo|'^^o-«o^o] ■ (37) 



dt2 

Finally, we have in the next order, O(e^), the averaged system 

z^— - Euo - a\uo\ Uq = - < K > vq + 2 [2\vo\ Uq - VqUq\ (38) 
zft— - Evq - a\vo\ vq^ - <K > vq + 2 [2\uo\ vq - UqVq\. (39) 

For periodic modulations of tunnel coupling we obtain M — Kf/Q'^. 

This averaged system is the main result of this section. We see that, in comparison with 
the constant coupling case, new effects appear. First, the cross term appears, corresponding 
to a change of type in the effective nonhnearity in coupled BEC dynamics. Second, there 
appears an additional nonlinear phase-sensitive coupling. This term leads to new minima 
in the effective potential picture. In the constant coupling case, the oscillations of z{t) is 
described by the Duffing oscillator. Now we have a more complicated nonlinear oscillator, 
describing more complicated interference patterns. 

V. ANALYSIS OF THE AVERAGED DYNAMICS 

Using the variables introduced in the section II, we get the system 

Zt^-2<K > sin($) - u5{l - z^) sin(2$) (40) 

$t = + 2 '^ ^ ^ ^ cos($) - u5z [2 - cos(2$)] . (41) 
Vl — 

where we have defined 5 = e^aM/h?. This system has the following Hamiltonian 

A 2 

H ^ -2< K > Vl- z^cos^- -{1- z"^) cos 2$, (42) 
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where Ai = A - 25. 

When the parameter 5 = 0, we get back the system (^. The effective particle motion is 
that of the double nonrigid pendulum system. 

Let us consider the different regimes of oscillations. 

A. Linear regime 

This is the case where 1^1 << 1,|$| << 1. Then, from the system (^OD, it follows that 
that 

= -2(< K > +5){K + 2 < K > -5)z. 

Then the frequency of linear oscillations of the number of atoms is 

ujI = -2{< K > +5){^ + 2 < K > -5) . 

Thus, the frequency of the linear oscillations grows with duo"^ = uj\ — uJq = 26{< K > +A — 5), 
in comparison with a nonmodulated case. 

B. z"^ << 1, and $ is not small. 

The phase dynamics is described by the equation: 

+ 2 < is: > (A - 26) sin $ + (5A + 2 < >^ -2(5^) sin(2<l>) + 

26<K> sin(3$) + — sin(4$) = 0. (43) 

At derivation of this equation we take into account that $f ~ z, so we neglect by terms ~ z^f 
This equation is equivalent to a generalized pendulum motion whose effective potential is 
given 

[/($) = -2<K>{A-25)-hA6 + 2<K>^ -25^) cos(2<l>) - ^ ^ cos(3$) - ^ cos 

2 3 8 

(44) 

The additional minima at $ = invr correspond to the trapping of oscillations in states witt 
locked phase. 
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VI. MACROSCOPIC QUANTUM SELF-TRAPPING. 



This phenomenon consists in the self trapping of the atomic number in one of conden- 
sates, depending on the initial difference in atomic number and on the relative phase of 
condensates . This is connected with the phase sensitive linear coupling. The nonlinear 
guided wave optics analog of this phenomenon, first predicted in , is the switching of the 
power in the dual core couplers from one core to other, depending on the initial power or the 
phase difference. The additional coupling in the averaged system (^) is nonlinear and also 
phase sensitive. Thus the condition of localization of atomic population should modified in 
comparison with the one considered earlier. 

Let us find the fixed points of the averaged system (^Oj). The first group is for = 
7m, n = 0, ±1, ±2, ... For each there can exist one or three roots for z. One of the roots 
is = 0, two others are 

^.--±(l-^X31^) . (45) 

Note that if (A — 5) < K >< 0, then 3 roots exists for even n and one root {zi = 0) 
exist for odd n. For (A — 5) < K » O,one has the opposite situation, the critical value is 
Ac> 2 < K > +5. 

One can show, however, linear stability analysis predicts that this set is unstable when 
6 > A/2. 

The second group of fixed points is given by: 

cos($) = -4^- (46) 



Vl-z 



Again, there are two possibilities: 

i) A = 36, then all line is stationary; 

ii) A 7^ 36, then one has the following static point 



cos($) = ^ = 0. (47) 
6 
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This root assumes that | < > /5| < 1. The hnear stabihty analysis shows that the fixed 
points are unstable when A > 36. 

Following , from the constraint coming from the energy conservation and the boundness 
of the tunneling energy, given by 

H = Ai^M -2«K» ^1-^2(0) cos ($(0)) - ^(1 - ^"(0)) cos(2<l>(0)) > 2 < K > 

(48) 

we obtain the estimate for the critical value of when the self-localization occurs. We 
find that: 



_ 4<K>{^1- z^jO) cos $(0) + 1) + ^((1 - z\0)) cos21>(0) - 1) 

The phase portrait is plotted in Fig. (9, 10, 11) for initial phase difference $(0) = vr and 
A = 2.5, for differ nt 6 = 0.2;0.6;0.8 The regions corresponding to different dynamics of z{t) 
with < z >= and < z >^ are shown. The increasing of S (i.e. Ki/Q) leads to the 
distortion of the MQST regime, and for 5 ~ 1 to the nonlinear Rabi-like oscillations. The 
typical oscillations in time of the relative population for different values of 6 and for fixed 
A are plotted in Fig.9. The infiuence of rapid modulations leads to the appearance of new 
minima in the oscillations. This can be described as a result of the overlap of two double 
well potentials. This conclusion is also confirmed by the results of numerical simulations of 



the GP equation for the single BEG under rapidly varying trap potential performed in |]18 
In this work, it has been shown that the effect of rapid modulations is the appearance of a 
effective double well trap potential. 



VII. CONCLUSION 

In this paper we have studied the new effects coming the time variation of the trap 
potential, and from damping, on the nonlinear oscillations in the relative population behavior 
between two BEG's. For the slowly varying trap we predict synchronization of oscillations 
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of the trap with oscillations of the relative population. We find the fixed points for both 
types of the damping terms occuring in the studies of two coupled condensates. Using the 
Melnikov approach we study the possibility the appearance of the chaotic oscillations in 
the tunneling phenomena between two coupled Bose-Einstein condensates. We calculate 
the width of the stochastic layer in phase space, in which the dynamics is chaotic. We 
find the lower bound on the region of chaos in {Ki,Q) plane. We obtain the estimate for 
the damping coefficient when the chaos is suppressed. For the rapidly varying trap we use 
the multiscale method and derive the averaged equations for coupled modes describing the 
tunneling phenomena. We find the fixed points in this case describing the stationary states 
in two coupled EEC's. The expression for the critical value of the nonlinearity parameter , 
when the macroscopic quantum self-trapping occurs is derived. Our results show that there 
is interest to study the contribution of the quantum tunneling for this driven case and the 
investigation of resonances and chaos in oscillations of atomic poulation for the strongly 



overlapped condensates p^. These problems require separate investigation. 
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FIGURES 
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FIG. 1. The typical behavior of the solutions to Eq.(^) when the damping is taken in the form 
= r](j)t. This case corresponds to = 0.7, r] = 0.1, Ki = 0.2, A = 2.5. Theoretical value is 
= 0.2599, the numerical is Zc = 0.26. 
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FIG. 2. The typical behavior of the solutions to Eq.(§) when the damping is taken in the form 
r = r]z. This case corresponds to = 0.7, AE = l,r] = 0.1, Ki = 0.2, A = 2.5 
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FIG. 3. The regions of regular and chaotic oscillations in plane Ki,Q from the Melnikov 
proach. 
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FIG. 4. The width of stochastic layer as a function of frequency when 
r] = 0,Ko = l,Ki = 0.2, A = 9.9, z{0) = 0.6, $(0) = 0. 
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FIG. 5. The typical oscillations of the relative population z{t) for 
Ki = 1.2, A = 9.9, z{0) = 0.6, $o = 0, = 1.5 
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FIG. 6. The phase portrait in the ($,2;)-plane for the same values of the parameters as in 
Fig(5). 
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FIG. 7. The typical oscillations of the relative population z{t) for 

= 0.2, A = 2.5, z{0) = 0.9, $o = tt, = 3.22 
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FIG. 10. The phase portrait in a plane z, (p of the averaged system of equations ( |40| ) for 
$(0) = 7r,A = 2.0,6 = 0.6 
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FIG. 12. The typical behavior of the solutions to Eq.@, for 6 = 0.8 and A = 10, $(0) = and 
varying z(0). (A) corresponds to z{0) = 0.5, (B) to z{0) = 0.65, (C) to z{0) = 0.68 and D to 
z{0) = 0.80. (B) is just before the localization transition and (C) just after. 
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